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Abstract. We calculate the bispectrum of single-field braneworld inflation, triggered 
by either an ordinary scalar field or a cosmological tachyon, by means of a gradient 
expansion of large-scale non-linear perturbations coupled to stochastic dynamics. The 
resulting effect is identical to that for single-field 4D standard inflation, the non- 
linearity parameter being proportional to the scalar spectral index in the limit of 
collapsing momentum. If the slow-roll approximation is assumed, braneworld and 
tachyon non-Gaussianities are subdominant with respect to the post-inflationary 
contribution. However, bulk physics may considerably strengthen the non-linear 
signatures. These features do not change significantly when considered in a non- 
commutative framework. 
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1. Introduction 

According to the inflationary paradigm, early-Universe small quantum fluctuations of 
the scalar field driving a period of accelerated expansion are amplified to cosmological 
scales by the expansion itself. These perturbations then leave their imprint in the 
cosmic microwave background (CMB) as thermal anisotropies. Two main physical 
observables are generated by this mechanism, namely, the scalar spectrum, which is 
(the Fourier transform of) the two-point correlation function of scalar perturbations, 
and the bispectrum, coming from the three-point function PP [Hj . For the curvature 
perturbation on comoving hypersurfaces TZ, this reads 

(^(k 1 )^(k 2 )^(k 3 )) = - (27T) 3 ^ (ki + k 2 + kg) 

x^2/ NL (k J ,k,)(|^ L (^)| 2 )(|^ L (^)| 2 ), (1) 

i<j 

where (|7?.l(^)| 2 ) is the power spectrum of the linear Gaussian part IZ^ of the curvature 
perturbation with comoving wavenumber k, sum indices run from 1 to 3, and /nl is 
the non-linearity parameter^. If /nl is momentum independent, one can write the 
gravitational potential in terms of 71^. in real space, 

ft(x) = ll h (x) - /nl [tt£(x) - <72£(x)>] , (2) 

which gives equation ((TJ) with / N l shifted outside the summation over the kjS. When the 
statistical distribution is Gaussian, /nl = 0, the three-point function vanishes. In terms 
of the CMB temperature fluctuation AT(e)/T, measured along the direction e, the limit 
of the bispectrum at zero angular separation is the skewness, Ss(e) = ((AT/T) 3 ). For 
practical purposes, this is a less sensitive probe for non-Gaussianity than the bispectrum 

It is common to restrict the discussion of cosmological perturbations to the power 
spectrum, implicitly assuming dealing with all the relevant informations that can be 
extrapolated from the sky. In many respects, the measured two-point correlation 
function is able, all by itself, to both describe the microwave sky in great detail and place 
observational constraints on the features of early-Universe models such as inflation. 

However, the experiments of the last generation have radically changed the general 
attitude towards cosmology and made possible what is now recognized as a 'precision 
era'. The physical scenarios explaining the large-scale structure of the Universe can be 
refined by more and more accurate observational inspections. Therefore it is natural 
to consider the bispectrum, too, and ask what signatures of non-Gaussianity we might 
expect from a given theoretical model (see [8 for a comprehensive review). 

A non-Gaussian spectrum can arise according to a number of different mechanisms: 
just to mention some, late-time non-linear evolution of cosmic structures [H] , multi-field 

X Our definition of the non-linearity parameter (sometimes dubbed /j^ L in literature) is —2 -3/5 that of 
(there denoted $3), —3/5 times that of and 3/5 times that of jElCj- The factor of 3/5 comes 
from the definition of /nl through the peculiar gravitational potential ^4, which is '5 4 = — 31Zl/5 
during matter domination. 
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inflation or inflation with scalar spectators ^U]-|I7|, the curvaton scenario [TH]-[20j. 
non- vacuum inflation |2I]~[2B]; higher-dimension operators in the inflaton Lagrangian 
[21], the inflaton with Dirac-Born-Infeld (DBI) Lagrangian |25jf. the ghost inflaton [26J, 
and, more commonly, a self-interacting inflaton |27| |31|. 

In four dimensions, it turns out that the inflationary contribution to the non-linear 
parameter is, for a single ordinary scalar field [H] E] , 



The scalar index can be written as a function of the slow-roll (SR) parameters describing 
the inflationary dynamics [321 E2] • At lowest SR order, that is in the quasi-de Sitter 
regime typical of inflationary models, the spectral index and the non-linearity parameter 
are first order in the SR parameters, which we denote symbolically as /nl = 0(e). 

It would be interesting to see how this result is modified for braneworld cosmologies 
and tachyon-driven inflation. In braneworld scenarios (e.g., [H3] and references therein) 
the visible universe and its matter content are confined into a (3+l)-dimensional 
Friedmann-Robertson-Walker (FRW) variety (a brane) embedded in a larger non- 
compact spacetime (called bulk). In the high-energy limit, the effective Friedmann 
equation experienced by a brane observer is, in suitable energy units, 



where H is the Hubble parameter, p is the energy density of the matter content (assumed 
to be a perfect fluid), and q is a constant. The Randall-Sundrum (RS) and Gauss- 
Bonnet (GB) braneworlds display a more complicated evolution, but in particular time 
intervals or energy 'patches' equation (0) is a good approximation, with q = 2 in the 
high-energy RS case and q = 2/3 in the high-energy GB case. Also, tachyon cosmology 
can be implemented in the same formalism with little extra effort. 

Intuitively, we do not expect a dramatic quantitative change in the non-Gaussian 
effect since the patch formalism does not interfere with SR expansions except for as 
regards the value of the coefficients in front of the SR parameters themselves |35| 
We have verified this guess explicitly. 

A preliminary confirmation comes from the stochastic approach of Gangui et al 
[3]. This approach permits one to estimate the order of magnitude of the effect by 
just considering second-order fluctuations of a self-interacting inflaton field and no 
gravitational fluctuations. This might seem too crude an approximation, since one 
should go up to second order in perturbation theory in order to fully take gravitational 
back-reaction into account and treat the bispectrum consistently. Surprisingly, the 
inflaton perturbation really encodes the main feature of the model apart from the 
resulting incorrect combination of SR parameters. 

f This model is quite distinct from the DBI cosmological tachyon. While the former predicts a strong 
non-Gaussian signature, the latter is more similar to standard inflation. 



4/nl ~ n s - 1 , 
where n s is the scalar spectral index 



(3) 




(4) 




(5) 
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Stochastic inflation is an approximated method according to which the solutions of 
the equation of motion for the scalar perturbation in the long wavelength limit k <C aH 
are connected to those in the k ^> aH limit at the Hubble horizon [37]— |42j. Let ip 
denote the scalar field on the brane. This (or its gauge- invariant version) is separated 
into a 'classical' or coarse-grained contribution ip( c \ encoding all the modes larger than 
the Hubble horizon, and a quantum or fine-grained part ip^ 9 ' taking into account the 
in-horizon modes. Therefore, the classical part is the average of the scalar field on a 
comoving volume with radius ~ (aH)' 1 . With this decomposition, the equation of 
motion for ijj becomes a Langevin equation with a stochastic noise source generated by 
the fine-grained contribution of the quantum fluctuations (see [13] for an extension to 
stochastic inflation with coloured noise). 

Actually, we have performed the calculation of jl] with the general FRW equation 
(JSJ) and the potential SR tower (SSj both for the scalar field and the tachyon. So, 
/nl — 0(e) for a generic braneworld filled with a scalar field. As we shall see, the 
key assumption giving rise to this behaviour is the validity of the standard continuity 
equation; see equation (fTlj) below. In fact, the presence of a brane-bulk flow would 
introduce genuinely novel features and considerably complicate the set of dynamical 
equations. Also, in the presence of extra dimensions the gravitational contribution 
may lead to a non-trivial behaviour of second-order perturbations, since to this order 
the interplay between extra-horizon scales and small scales may become quite delicate. 
This would impose a more rigorous treatment and a full second-order calculation in 
order to carefully evaluating the non-Gaussianity produced during inflation, e.g., one 
like that performed in jHj for the four-dimensional case. 

A powerful approach in this sense is the space-gradient formalism of jS]-|17J, a 
development of the separate universe method 0H] (see also the earlier works |4^]-|53j). 
We shall follow the notation of the cited papers closely, and skip details that can be 
found there unaltered by brane or tachyon physics. 

In this paper we achieve three goals. The first is to demonstrate the validity of 
the Mukhanov equation for the braneworld at large scales (section EJ); this equation 
allows us to calculate the power spectrum of the gauge-invariant scalar field fluctuation 
u oc a5ip through a SR expansion [Mj-jEE], and was long suspected to be valid even in 
the braneworld case. Not only does the method that we shall adopt give independent 
support to other proofs (e.g., jSTJ-jSH] for the RS brane), but also it is valid beyond 
the linear order in perturbation theory. Also, our equation holds for a general patch 
and not only for the Randall-Sundrum case. The second outcome is the bispectrum 
of perturbations generated by either a cosmological tachyon or scenarios, such as high- 
energy braneworlds, with a modified Friedmann equation (section HJ); so far, only in 
|35| lol?] have some (very preliminary) remarks been made on the subject. Finally, 
in section 03 these results are extended to the case of a maximally symmetric non- 
commutative spacetime jHO]- Non-commutative cosmology has recently drawn a lot of 
attention due to its possible effects on the shape of the power spectrum (for a list of 
references, see jHUE^l); here we shall inspect for the first time the bispectrum predicted 
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by this class of models. Conclusions are given in section El 



2. Set-up 

To begin with, it is convenient to work in the induced metric 

ds 2 4 = N 2 (t, x)dt 2 - a 2 (t, x)dx i dx i , (6) 

brane 

where N(t,x) is the lapse function (TV = 1 in synchronous gauge), a(t,x) is a locally 
defined scale factor, and {x%, X2, x^} are space comoving coordinates on the brane. The 
Hubble parameter reads H = a/(Na); dots (primes) will denote derivatives with respect 
to t (ip). Also, we define II = ip/N and the parameter 9 = 2(1 — q^ 1 ). The first three 
SR parameters, which will be all we need, are 

~ NH 2 2 H 2 ' { ' 

ri 

(8) 



NHTL 



^ ~ nh 2 (/vn J ' (9) 

where 9 = 9 for the ordinary scalar and 9 = 2 for the tachyon. The SR evolution 
equations are 



NHe 



277 



(10) 

f} = NH(er]-t 2 ) . (11) 

In the slow-roll approximation, e, |r/|,^ 2 <C 1, the inflaton potential energy V(ip) 

dominates other the kinetic term and the accelerated expansion is quasi-de Sitter 

(H,e,rj ~ constant). Sometimes a subscript in the SR symbols will specify the scalar 
model. For an ordinary scalar field 0, the equation of motion at large scales reads 
(IT = 0//V) 

il + 3NHTI + NV'((f)) = , (12) 
while for the DBI tachyon T one has (II = T/N) 

U tt0 + 3NH11 + N\nV'(T) = . (13) 
1 — IT* 

These are nothing but the continuity equation 

p + 3NHp(l + w) = 0, (14) 

where the barotropic index w(p, p) relates the pressure and energy density of the perfect 
fluid filling the brane, p = wp. 

In the separate universe approach, the physical quantities such as H(t,x), the 
scalar field ^(t, x), cosmological perturbations and observables are defined on an 
inhomogeneous background and evolve separately through the dynamical equations at 
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each point once the initial conditions have been specified. Then we can transform time 
derivatives into spatial gradients like 



h n r ' 



0~l)e- V 



§<9^, (15) 



where the index i runs from 1 to 3, 9 = 1 for the ordinary scalar field, and 6 = 8/2 for 
the tachyon. 

At large scales, second-order gradient terms can be neglected and the gradient 
generalization of the Bardeen potential (HHI is conserved |44j : 

C . = Xi -NH^, (16) 
P 

6 = 0, (17) 

where Xi = <9jlna. In the braneworld case, the long wavelength limit is advocated for 
consistently neglecting the projected Weyl tensor. This in turn is deeply intertwined 
with the other fundamental constraint, that is a standard continuity equation (|14j) . In 
the presence of a source term surviving at all scales in the continuity equation, equation 
(fT7jl would not be valid and the quantity Q would not be conserved. However, bulk 
physics affects only the small-scale or late-time cosmological structure and does not 
play a chief role during inflation. This result has been confirmed with several methods 
both analytically and numerically [IB], |54"j-|68j. We shall come back to this issue in 
the final section^. 

For a scalar field, the curvature perturbation ( on hypersurfaces with constant 
energy density coincides with the curvature perturbation TZ on comoving hypersurfaces, 
which as a vector quantity reads 

Ki=Xi- ^drf . (18) 
In the linear approximation, equation (fTH|) is the spatial gradient of 

where \l/ 4 = —5a/ a is the gauge-invariant potential perturbing the spatial part of the 
metric. Define 

an / 2e \ 1/2 , , 

* s ^ = • Up*) • < 20) 

I Also, the evolution equations for the SR parameters defined by equations J7J and ©, as well as 
all the expressions which will be derived through equations lfTT)|l and JTJJ, depend on the continuity 
equation with no extra flux. However, one can define the SR parameters only in terms of the energy 
density with no reference to equation l|14[) . Let N = L H(t')dt' be the number of e- foldings counted 
forward in time, and a subscript TV denote differentiation with respect to it. Then, for instance, the 
parameters e(p) = —qpx/(2p) and rj(p) = —pNM '/(2pAr) satisfy equation (| 1 for an ordinary scalar 
field, and coincide with the usual SR quantities only when equation (|14|l holds. 
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where the speed of sound cs is cs = 1 for the ordinary scalar field and cs = 
v/T— 2er/(3g) for the tachyon. The spatial gradient of the Mukhanov-Sasaki variable 
u = —zTZl corresponds to the linear limit of 

Qi = -zd = -zKi = a (d4 - ^X^j , (21) 

where a = zH/U. For the ordinary scalar field, a = a. From equations (jl5j) . ()20|) and 
(ED) we get 

diZ = -Qi + ad^[l + 9e-r] + (e- l)^ 2 ] , (22) 
where the quantity 

(23) 



2 _ q c s _ 4ge T ?7 T 



NH c s 3g - 2e T 
is second order in the SR parameters. 

3. Generalized Mukhanov equation and stochastic inflation 

The evolution equation for can be computed within the multi-field framework of 
[4*H] . Here we consider the same calculation for a single field ip in a patch given by 
equation (J3J); the results will match each other in the 4D case. The first time derivative 
of equation (f2~Tj) is 

Q l = - z Qi = NHQi [l + 0e-ri + (6- l)^ 2 ] . (24) 
Another derivation gives 

Qi--Qi = 0. (25) 
However, it is more convenient to recast the equation of motion as 

Q.-FQi + nQ, = 0, (26) 

where 

F^^-NH, (27) 

A = p- - - = 2 + (30 - 1) e - 3?7 - 4^ + (l + - PjOe 2 + r? 2 + £ 2 + (6 - l)g^ 2 - 

(28) 

The extra tachyonic term is 

The expression for f2 is in accordance with the computation in the linear theory jSS] and, 
as that, is exact in the SR parameters. The expression found in jJHJ EE] (0 = 0,6 = 1) 



Qq = 3 + (Q _ i) eT _ 2r/T 4 ( - - 1 ) ^ 2 + ^2" • ( 29 ) 



is recovered via the exact relation 

,2 (-2 



V" 



H 2 



3(e* + 7fe)-74-£j. (30) 
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Equation (J26)) is equivalent to a generalized Mukhanov equation (with the Laplacian 
term dropped) when expressed via conformal time dr] = Ndt/a; since d 2 = (N/a) 2 d 2 + 
(N/a)Fd v and = — (N/a) 2 d 2 z/z, one has 

d 2 ld 2 z\ 

drj 2 z drf ) ^ ' 

In the linear approximation and in momentum space, equations (j26|) and (}3T| hold 
for Qi pa diU — > i^Mk, where a subscript denotes the Fourier mode with comoving 
wavenumber k. A standard calculation in de Sitter spacetime gives the mean value of 
the quantum field w k : 

aH 



(32) 

The equation of motion for can be written as an equation for the coarse-grained 
part of «k sourced via a stochastic noise term. The coarse-grained part of the Mukhanov 
variable is u k = u\ s yV(kR), where W is the Fourier transform of a window function 
W(|x — x!\/R) falling off at distances larger than R. The scale R is of order of the 
comoving Hubble radius, R = /i(aif) -1 , with h > 1 so as to encompass the whole 
horizon. With h sufficiently larger than 1 we can safely discard the k 2 term in the 
Mukhanov equation. 

The final result is then extended to the non-linear gradient variable Qi at large 
scales, yielding 

Qi - FQi + flQi = £i(t, x) , (33) 
£(t, x) = / ^^i^e ik - x W + W(2d t - F) n k a(k) + c.c. , (34) 



(2tt) 3 /2 

where the superscript (c) in Qi is understood and a(k) is a complex stochastic quantity 
such that the ensemble average (a(k)a*(k')) = 5^ 3 ^(k — k'); it simulates the continuous 
crossing of modes outside the horizon and their addition to the coarse-grained part. As 
they stand, equations (J3~3*j) and properly encode the full stochastic contribution. 

Equation (jHSJ) is the non-linear extension of the Langevin-type equation we used in 
the first- SR-order heuristic computation in synchronous gauge: 

^ , w-w-l5 + « x ' t,ss -i? +5(x ' t) ' (35) 

where U(<f>) = V(<p) and U(T) = In V(T). In this and the other equation one can see 
that there are basically two sources of non-linearity. The first one is the back-reaction 
of the field fluctuations on the background encoded in the noise term; the second one 
is the self-interaction of the scalar field represented by the potential contribution Q 
(or —U'/(3H)). Therefore, a priori the statistical distribution of Q^ (-0^) will be 
non-Gaussian, even if quantum fluctuations are random. 
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4. Braneworld and tachyon bispectrum 

In order to compute the scalar spectrum and bispectrum, we fix the gauge to the time 
slicing with respect to which the kth mode crosses the horizon simultaneously for all 
spatial points jH]. Then t = ln(aif), NH = (1 — e) -1 , and R = he~ f . In this gauge, 
the gradient curvature perturbation and diZ are, respectively, 

Qi = ad^{\ - e) , (36) 
3z« [(l + 9)e-r{\ Qi (37) 
= -l{n s -l)Qi, (38) 

where in the last passage we have used the first-order patch SR expression for the scalar 
spectral index |331 l3~o]: 

n s -l = 2r}- 2(1 + 6) e. (39) 

At first SR order, 

F ~ —1 , 2 + 3(9 + 1) e — 3»7 . (40) 

Equation can be expressed as a Langevin differential equation in the curvature 
perturbation: 

Hi + (2 Z - - F^j TZi = ~ £i(t, x) . (41) 

To lowest SR order and neglecting the IZi term, 

^w-^ftx), (42) 

where cs ~ 1 inside & and z. The power spectrum is given by the solution of the 
linearized equation (|42|) : at first order in a perturbative expansion, 

-m 1 ' « — - (43) 

where z^ is z defined on the homogeneous background. The time integration of £/<z from 
the initial time t t to t is proportional to B(kR) = j£ d^a#) _1 [>V + VV(2d t , + l)](a#) = 
[1 + (kR) 2 /3)W{kR) - [1 + (kRif /3}W{kRi) to lowest SR order. Here we have used a 
Gaussian window function 

W = exp(-/e 2 i? 2 /2) , VV = (A;i?) 2 W , W = [(A;i?) 2 - 2}W , (44) 

and the lowest- SR-order eigenvalue equation d t Uk = u^. In the limit of asymptotic 
past (ti — > — oo) and future (t — > +oo), B(kR) — > 1. The integration over fc yields 
= T^ 1 ^ = d~ 2 d l lZi and the lowest-order amplitude 

<|K<1,w|2>=3 ifV- (45) 
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after a computation almost identical to that of [IE]. At second order in the perturbation, 
-mf ] = (z- 1 )^^ + 4 {2) /^ (0) - Since ^ (2) = 0(e 2 ), the only surviving term at first SR 
order is, by the first-order version of equation (JHHJ), 

= -l(n B -l)n^\ (47) 

After computing the non- linear term ft (2) = -±(n s - 1)[^ (1) ] 2 , one is ready to write 
down the curvature perturbation at second order, 

K « + ±ft (2) , (48) 

and the bispectrum Actually, the integration of the window function inside the noise 
term corresponds to a non-linearity parameter with non-trivial momentum dependence, 
for which equation (J2J) does not hold. When one of the momenta is negligible relative 
to the others, that is in the squeezed limit k 3 <C ki,k 2 , by equations (j4Tj) and (JTJ) we get 

4/nl = K - 1) + 4/(k 1; k 2 , kg) « n s - 1 . (49) 

Although we have not written explicitly the momentum structure f(k), which at first 
SR order is the same as in [46J, we can draw some important conclusions. 



(i) Tachyon and ordinary inflation generate the same non-Gaussian signature in 
the limit of collapsing momentum dependence and at first SR order. Outside 
this approximation the non-linearity parameter /nl can depend on the type 
of scalar field and braneworld through the parameter z: 4/nl = { n s — 1) + 
4/( fl '^(e, 77; ki, k2, kg). This is in agreement with the correspondence between the 
lowest-SR-order tachyon and ordinary observables. 

(ii) Written in terms of n s , braneworld non-Gaussianity does not differ from the 
4D picture, except perhaps in higher-order contributions. What changes is the 
inflationary model one has to impose in order to predict a given scalar spectral 
index. In this sense, we could regard equation Q as a first-SR-order consistency 
equation joining the traditional set we explored elsewhere [HE], since both n s and 
/nl (through the bispectrum) are observables [701 ITT] $. 

(iii) On the other hand, one should note that the post-inflationary era greatly enhances 
non-Gaussianity, up to /j^ = 0(1) [Z2]^[Z11|| or even /^L St ~ 50 from a 
suitable preheating phase [ZHj. As explained in |Hj, the true observed non-linearity 

§ A small non-Gaussian component comes also from the three-point functions involving the graviton 
zero mode. Using the z function for braneworld tensor perturbations z cx aH~ 8 / 2 |8f>| , one finds a 
contribution proportional to the tensor amplitude and spectral index n t — — (2 + 0)e, in agreement 
with the 4D result jxj- However, since the tensor amplitude is much smaller than the scalar one, we 
can neglect this term with respect to the scalar bispectrum. 

|| This is basically due to the fact that at second order the longitudinal gauge condition $^ — = is 
(2) (2) ( {\\\ 

modified as <9\' — = 4 I ty y 4 ' \ at large scales, thus providing a non-trivial second-order correction 
to the Sachs- Wolfe effect 0- 
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parameter is not the bare inflationary result (J3J). In addition to the post-inflationary 
contribution, one must consider angular averaging. The total observed /nl is in 
fact, and at least, = 0(1)+/nl- Therefore the non-linear effect of braneworld or 
4D inflation, if the SR approximation holds as we required, is always subdominant. 

Although these features might be obvious when inspecting patch cosmology, here we 
have derived them quantitatively^. Moreover, the gradient+stochastic approach can 
be the basis for next-to-leading-order SR and perturbation calculations as well as for 
numerical simulations [35] . 

According to the first-year WMAP analysis, the power spectrum fully characterizes 
the statistical properties of CMB anisotropies: that is, /nl vanishes consistently. More 
precisely, the constraint on the non-linearity parameter is —58 < /nl < 134 [76J, 
which does not discard not only inflationary non-Gaussianity, but also other models 
predicting a more robust effect, / NL 3> 1. See also [77J[7H] for other analyses. The next- 
year WMAP data and the Planck satellite should significantly improve the accuracy of 
the measure, with the inclusion of polarization anisotropies: /^"(WMAP) ~ 11 — 15, 
/^(Planck) ~ 3-5 [13]. 

We also note that in this class of models (single brane field and empty bulk) the 
contribution of f(k) is below the present observational threshold. Actually, in order to 
experimentally detect and constrain the contribution of the momentum-dependent part 
°f /nl, h should be much larger than 0(e), say f(k) ~ 0(10). Since the characteristic 
coefficient in f(k) is 0(1) (see [IB] for details), the exact form of this term is not relevant 
in such a set-up. 

5. Non-commutative bispectrum 

Until now we have considered a commutative background throughout the whole 
spacetime. We can make a step further and phenomenologically assume that we have a 
3-brane in which the stringy spacetime uncertainty relation (SSUR), At Ax > 1% , holds 
for all the braneworld coordinates {x u }, v = 0, 1 2 3, while the extra dimension y along 
the bulk remains decoupled from the associated *-algebra. Here, r = J dtNa, (« a/H 
in the SR regime), x is a comoving spatial coordinate on the brane, and l s = M~ l is the 
fundamental string scale. 

Non- commutative braneworld inflation arises when we impose a realization of the *- 
algebra on the brane coordinates. In order to diagonalize the non-commutative algebra 
and induce a pure 4D SSUR on the brane, one might fix the expectation values of the 
background fields of the fundamental theory such that the extra direction commutes, 
[t/, x^] = 0. An algebra preserving the maximal symmetry of the FRW universe is 
[t,x] = i/g [60J, by which normal products in the action of the perturbation modes are 

% The authors of |71| drew similar conclusions, claiming that equation Q is a model-independent 
consistency equation under the assumption of single-field inflation and in the squeezed limit. However, 
see [H] (Sec. 8.4.1). Reference [23] deals with the de Sitter case only. 
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replaced by *-products like 



(50) 



The non-commutative algebra induces a cut-off k roughly dividing the space of 
comoving wavenumbers into two regions, one encoding the ultra-violet (UV), small-scale 
perturbations generated inside the Hubble horizon (H <C M s ) and the other describing 
the infra-red (IR), large-scale perturbations created outside the horizon (H 3> M B ). By 
definition, they correspond to the quasi- commutative and strongly non- commutative 
regime, respectively. 

The separate universe approach does not contrast with a non-commutative 
background. To understand this point, we can use the linear picture of |48j . and in 
particular their figure 1. The basic idea is that a comoving large-scale perturbation 
is independently specified in two comoving locally homogeneous regions separated by a 
distance A, if the size A s > H~ l of these regions is small with respect to A. Perturbations 
are defined on a given background, that is a region of scale Ao much larger than our 
present horizon. Then the required hierarchy of scales is Ao ^ A ^> A s > H -1 . In the 
presence of a non-local algebra, the string scale can play the role of natural marker in 
the hierarchy. For instance, setting A s ~ l s we just consider the IR region of *-models. 
If A s > H^ 1 > l s , the previous argument is unchanged. 

Since the *-product (|5U|) does not involve homogeneous quantities (i.e., it preserves 
the FRW maximal symmetry), the Mukhanov equation for a non- commutative 4D (HHj 
or braneworld [SI] scenario is, at linear order and large scales, 



where the effective conformal time drj = (a/a e s) 2 dr], effective scale factor a e g, and 
measure Zk are constructed through smeared, momentum-dependent versions of the 
scale factor a± = a{r ± kl 2 ). In the separate universe approach, a± acquires a spatial 
dependence like the other quantities, a±(r) — > a±(x, r). + The measure z k is given by the 
product of z and a correction factor f z depending on the particular non- commutative 
model one is assuming |H] . The de Sitter solution of equation (JHT|) is the commutative 
solution (JH2*j) multiplied by fl = (a c{i /a) 2 , which is the relative rescaling of commutative 
to non-commutative conformal time. 

The above discussion on conserved non-linear perturbations is not modified by 
the introduction of a fundamental length scale. Therefore we are tempted to directly 
generalize equation (f5*Tf with the gradient variable Q,. Since the mass term d~Zk/zk 
now depends on k, we define the Fourier mode of Qi as 



+ Note that the correct procedure is first to smear the scale factor in a sufficiently small homogeneous 
neighbourhood, a(r) — > a±(r), and then to extend it to very large scales, a±(r) — > a±(x, r). The 
'top-down' smearing a(x, r) — > a±(x, t) does not lead to equation (|5f (I . 




(51) 




3 



(52) 
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Then equation (joT|) holds for w k — > <2(k). One gets the t- version equation (pjj) by using 
the effective time dt = drja c g/N = dt/f a and setting N = N without loss of generality: 

d|Q(k) - FdfQ(k) + fiQ(k) = , (53) 

where F and f2 are 

F = f a F-f a , (54) 

n = - ^ = - f — ) 2 ^ = - u,) , (55) 





(56) 



In the commutative limit, F — > F, u — > 0, and — > Q. The term in square brackets in 
equation (jSUjl contains the 0(e) contribution of fl, since / Q // a = 0(e) = f z /fz- In the IR 
limit (strongly non-commutative regime), to and its components loose their momentum 
dependence [fJT]. The Langevin equation for the curvature perturbation reads 

dffc(k) + ( 2f a ^ - dfW(k) = -- k) . (57) 

V z k J z k 



The non-commutative version of equation (|22jl gets an extra term from the redefinition of 
z; in momentum space and at first SR order, Zk ~ — Q(k) + akipk[l + 8e — rj + f z /(NHf z )}. 

With the gauge choice rj pa — (aif/ 2 ) -1 = — e - *, by definition one has d{t = on 
surfaces of constant non-commutative time. Then 

NHf a = = i + 0(e) , F = -1 + 0(e 2 ) , d^u k = u k . (58) 

1 — e 

Neglecting the second-derivative term in equation 1)5 7 j). one has 

dfW(k) «--?-£(*, k). (59) 
3z fc 

With the procedure of the last section, at lowest SR order and first perturbative 
order, one obtains the non-commutative power spectrum found in |6TH I5T] . which is 
the commutative amplitude corrected by a factor £ 2 = (f%/ f z ) 2 - Actually, if one 
wanted to go to coordinate space, the integration over momenta should be performed 
up to the UV cut-off fc , while the characteristic scale R should be pushed towards the 
asymptotic limit R ~ k^ 1 at the infinite future. However, the approximation ko — > oo 
fits well at this stage of accuracy. 

The gauge-fixed variables Q and Zk are 

Q(k) =a k ^ k (l-e + 2/ a ), (60) 



Zk 



{l + e)e-r 1 + fj-f-2f a 



Q(k) = -^K-l)Q(k). (61) 



At first SR order, the spectral index ()39|) has acquired an extra term 

dlnS 2 _ 2 
dlnfc ~ ~NH 





ae = = -A- i 2^ - ^ ! « 2,/;, I 2^ - ±L I . (G2) 
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where a is constant in the IR limit. A computation of the second-order curvature 
perturbation, by equation (JoTj) . gives again equation Q, with the spectral index now 
depending on the non-commutative parameter a. 

6. Conclusions 

In the context of tachyon and braneworld scenarios, early-Universe large-scale 
perturbations can leave an imprint different to that of standard four-dimensional 
inflation in the microwave sky. Under certain assumptions, we have shown that the 
linear cosmological spectrum comes from the first term of a gradient perturbative 
expansion of a non-linear curvature perturbation satisfying a generalized Mukhanov 
equation of motion. The bispectrum of this quantity, which involves it at second order, 
governs the non-Gaussian signature eventually detectable in the CMB. Neglecting the 
projected Weyl tensor on the brane, we have found that the non-linearity parameter /nl 
is proportional to the scalar spectral index and therefore unobservable, in agreement 
with past 4D calculations and even in the case of a particular realization of a non- 
commutative spacetime. 

A way to generate a stronger non- Gaussian signal in the braneworld context might 
be to include the Weyl contribution (both in the Friedmann and continuity equations), 
but as we have said we should expect it to play a negligible role in the long wavelength 
limit. Nevertheless, this issue will deserve further attention for at least two good reasons. 
The first is that Weyl damping was considered and shown only in the case of linear 
perturbations, while the stochastic Langevin equation for the curvature invariant holds 
at all orders. The second is that bulk physics intrinsically provides a noise source to the 
Mukhanov equation through an infinite tower of Kaluza-Klein scalar modes dominating 
at short wavelengths |HB]- Therefore, while equation would keep on being valid, 
an important contribution to the stochastic noise term (JMj) might be lacking in the 
present analysis. In this case the consistency equations, included equation (JHJ) even in 
the squeezed limit, would be spoiled anyway 

Therefore one should expect to find the most interesting non-Gaussian effects, 
if any, in the spectral regions we have neglected so far, namely, around the time of 
horizon-crossing or after, at sub-horizon scales. Here one should note a caveat hidden 
in the above approach. The stochastic noise given by the random field a(k) mimics at- 
horizon and sub-horizon physics without a more fundamental knowledge of it, so that 
this method implicitly misses cubic interactions around horizon crossing. As stressed 
in jHU], this might be the source of the discrepancy between jlH] and [7] as regards the 
shape of f{k). In the latter paper this physics is under control and the resulting answer 
perhaps relies upon more solid ground. For these reasons, one might be concerned with 
whether the formalism we chose is not the most suitable to investigate signatures of 
braneworld non-Gaussianity. However, it is probably one of the most calculationally 
accessible methods at present time, which gives the correct momentum-independent 
result anyway. In a more complete treatment of braneworld non-Gaussianity, the k- 
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dependent part of the non-linear parameter might change relative to the 4D result and 
should be presented explicitly, at least for estimating its magnitude. As stated in section 
HJ only when f(k) > 10 could this term be important from the observational point of 
view. Only if this were the case in the stochastic approach would the search for an 
alternative method be worth the extra effort. 

To conclude, we would like to be more tentative in guessing the magnitude of 
/nl when the bulk physics is taken into account. In general, the effect is rather model 
dependent and there is still no univocal result. For example, large non-Gaussianities can 
be obtained in the 'D-cceleration' scenario [23] or for specific choices of the Weyl tensor 
in the RS brane (HS] • However, in [HI] the Weyl noise term in the M ukhanov equation for 
the RS brane has been calculated and is 0(e) (see their equations (103)-(106)), which 
might imply a small extra contribution to the three-point correlation function, not larger 
than the spectral index. This issue should be checked and is now under investigation. 
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